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In this paper we consider cases of existence of invariant measure, additional first integrals, and Poisson structure in a 
problem of rigid body's rolling without sliding on plane and sphere. The problem of rigid body's motion on plane was 
studied by S. A. Chaplygin, P. Appel, D.Korteweg. They showed that the equations of motion are reduced to a second- 
order linear differential equation in the case when the surface of dynamically symmetric body is a surface of revolution. 
These results were partially generalized by P. Woronetz, who studied the motion of body of revolution and the motion 
of round disk with sharp edge on the surface of sphere. In both cases the systems are Euler - Jacobi integrable and have 
additional integrals and invariant measure. It turns out that after some change of time defined by reducing multiplier, 
the reduced system is a Hamiltonian system. Here we consider different cases when the integrals and invariant measure 
can be presented as finite algebraic expressions. 

We also consider the generalized problem of rolling of dynamically nonsymmetric Chaplygin ball. The results of 
studies are presented as tables that describe the hierarchy of existence of various tensor invariants: invariant measure, 
integrals, and Poisson structure in the considered problems. 



> 
•i—t 

X 



Contents 



1. Equations of Rigid Body Motion on Plane and Sphere without Sliding (Nonholonomic Rolling}] ! 1 



2. Body on a Plane! 180 

Body of revolution's rolling on plane (S. A. Chaplygin P. Appell [1, 131). 180 



Round disk (S. A. Chaplygin. P. Appcll. D. Korteweg).! 183 



Dynamically symmetric ball with the displaced center of mass (E. Routh. S. A. Chaplygin) 



184 



Three-dimensional point maps in nonholonomic mechanics.1 184 



Rolling of balanced, dynamically nonsymmetric ball (Chaplygin ball 125 ). 1 185 



Rolling of unbalanced, dynamically nonsymmetric ball on plane J 186 

An arbitrary body with a spherical central ellipsoid of inertia] 187 



Gvrostatic generalizations 



188 

R olling of ellipsoid on plane J 189 



Mathematics Subject Classification 37J60, 37J35 



REGULAR AND CHAOTIC DYNAMICS, V. 7, A*2, 2002 



177 



A. V. BORISOV, I. S. MAMAEV 



3. Body on a Sphere! . . . 190 

Rolling of body of revolution 192 



The problem of rolling of round disk! 193 



Ball with displaced center 



193 



Balanced, dynamically nonsvmmetric ball on sphere.1 . 194 



Unbalanced, dynamically nonsvmmetric ball on sphere.! 194 



Gvrostatic generalizations] 194 



Rolling of body with partially flat surface on sphere 



195 



Conclusions! 197 



1. Equations of Rigid Body Motion on Plane and Sphere without 
Sliding (Nonholonomic Rolling) 

In this paper we consider the equations of rigid body's rolling on 
plane and sphere, because in these two cases, as opposed to a 
rolling on an arbitrary surface, the equations of motion are sim- 
ilar the Euler - Poisson equations. In both cases there are six 
first-order equations for six variables. In the potential field the 
equations have two integrals of motion: the energy integral and 
the geometrical integral (for Euler - Poisson equations there is also 
the area integral; its analog in problem of rolling is presented only 
under additional dynamical and geometrical restrictions). 

Suppose that the rigid body rolls without sliding (i.e. the 
velocity of contact point Q is equal to zero) on the fixed surface 
represented by plane or sphere. The first part of equations of 
motion is the vector dynamical equation of kinetic moment M 
behavior in time with respect to the contact point Q (Fig. This equation is represented for 
arbitrary shapes of body and surface in the form 




Fig. 1 



M = M x u + mi- x (u x r) + Mr 



(1.1) 



where M, u, r = GQ, Mq are supposed to be projected on the principal central axes of inertia in 
the body; here u> is the angular velocity, Mq is the moment of external forces with respect to the 
contact point, G is the center of mass. The second part of the motion equation is the vector kinetic 
equation of Poisson type different for the cases of plane a) and sphere b): 

a) 7 = 7 x w, (1-2) 
where 7 is the unit vector orthogonal to the plane, 

b) R (-y + u x~f) = r, (1.3) 
where 7 is the unit vector orthogonal to the sphere of radius Rq (see Fig. |HJ). 

In equations (|1.1|) , (|1.2|) , Q1.3JI we suppose that the radius vector r is expressed as a function of 
normal vector 7 with the help of equation 



7 



grad/ 
grad/| 



(1.4) 



that define the Gauss transformation, where f(r) = is the equation of the body's surface in the main 
central frame of references connected with the body. We suppose that the body is everywhere convex 
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(to exclude the collisions during the motion), and equation (J1.4JI is uniquely solvable with respect 
to r = r('y). We assume that this condition is always fulfilled in the following text. 

Using (jl.4j) we can represent kinematic equation (jl.3j) describing the dynamics of vector 7 in the 
case of sphere in the form 

7 = (1 + fc(B - k)- 1 ) 7 x u, k = R~\ (1.5) 

where B = ||6,-,-|| is a degenerate matrix with the components ba = — { — 

11 u " 6 L 13 dn VjV/| drj 

The relation between M and u is defined by the equation 

AM — mr x (Ar x AM) 

M = Iti> + mr X (u> X r), cj = — , (1.6) 

1 — m(r, Ar) 

where m is the mass of body, I is the central tensor of inertia, A = (I + mr 2 ) -1 . 

If the potential U = U (7) depends only on the components of vector 7, then we can present the 
moment of external forces in the form: 

or r 

a) in the case of plane Mq = 7 x ; 

b) in the case of sphere Mq = 7 x (1 + k(B T — k)' 1 )^-, where B is matrix l|1.5j) . 

For equations (|1.1|) and (|1.2|) , Q1.3JI we always have the energy integral and the geometrical integral 

£T=i(M,w) + tf(7), F 1 = 7 2 = l. (1.7) 
Remark . The proof is based on the formula 

i(Af > w)- = (M ,w). 

It follows only from equation (f 1 . L|t and does not depend on the shape of the surface, on which the body is 
rolling. 

According to Euler-Jacobi theorem (the theory of last multiplier), to integrate these equations 
we need two more independent first integrals and an invariant measure JS]. Recall that the density 
of invariant measure p of the general system 

x = v(x), x = (xi, . . . ,x n ) 

satisfies the Liouville equation div(/?u) = 0. In the general case none of these objects exists; therefore, 
the system shows some interesting asymptotic and chaotic properties characteristic for the oscillation 
of Celtic stones. Such properties are not typical for conservative systems [131 115j . We consider all the 
known cases of existence of additional first integrals (one or two at once) and the cases of existence of 
invariant measure. 

Dynamical and geometrical bounds leading to the existence of the first integrals and of the 
invariant measure are independent in some sense. For some combinations of parameters there exist 
only the measure or only the additional integral. In some extreme cases two additional integrals and 
measure can exist at the same time; thus, the system became completely integrable. 

The results of the study are presented separately for the cases of body's rolling on plane and on 
sphere and collected in tables 1, 2. The following subsections are essentially the comments for these 
tables. 
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2. Body on a Plane 

Body of revolution's rolling on plane (S. A. Chaplygin |26j, P. Appell [1, 31J). If both 
the surface of body and the central ellipsoid of inertia are coaxial surfaces of revolution, then for 
equations 1)1.1(1 . 1)1.2)1 there exist two additional integrals and invariant measure. We assume that the 
potential U is an arbitrary function of 73 = cos#, i. e. it depends only on the slope of revolution axis 
of the body to the vertical. In particular, in the case of the gravity field the center of mass must be 
situated on the axis of revolution. 

The integrability of problem for an arbitrary body of revolution was shown by S. A. Chaplygin in 
1897 [2n]- He also demonstrated that it is possible to add a balanced uniformly revolving rotor along 
the axis of revolution (gyrostat) preserving the integrability of problem. More specific cases of this 
problem were studied by 

a. Routh (1884): the rolling of unbalanced dynamically symmetric ball on plane. 

b. Neumann, Carvallo (1898), Appell (1899) and Korteweg (1900): the rolling of round disk. 
The results of Neumann and Carvallo mainly concern the deduction of motion equations and 

determination of stationary solutions. Note that Neumann during the deduction of equations of motion 
at first made the same mistake that occurred before in Lindelof paper. He applied the Lagrangean 
formalism without the necessary " nonholonomic" modifications. In the subsequent studies he corrected 
this mistake, but did not solve the problem in quadratures. The lindelof mistake was analyzed in detail 
by S. A. Chaplygin (1897). He obtained a new form of equations of nonholonomic dynamics and was 
able to reduce the considered problem of rolling of revolution body to two linear first-order equations. 
In the case of round disk's rolling, S. A. Chaplygin showed the possibility of reduction of these two 
equations to one linear (second-order) equation solvable in hypergeometric functions. We should also 
note that, before Chaplygin's work, the equations of heavy revolution body motion were obtained 
in 1861 by G. Slesser, but their integrability was not indicated. 

Somewhat later (in 1898), the analogous substitution (in the equations obtained by Carvallo in 
the paper presented for the Fourneyron prize) was used by Appell and in the slightly different form by 
Korteweg. They both did not know S. A. Chaplygin's paper that was published in inaccessible journal 
only in Russian (the English translation of Chaplygin's papers dedicated to this problem is published 
in this journal in 2002). This is the reason of the fact that in many modern textbooks and research 
papers (O'Reily jSHj) the problem of round disk's rolling is connected with the names of Appell and 
Korteweg, although the previous text show that this opinion is not completely correct. 

Here we present the results obtained by S. A. Chaplygin in the modern algebraic form that let us 
to show the invariant measure in the explicit form and also to obtain the simplest forms of the first 
integrals. It turns out that we can generalize these results to dynamically nonsymmetric situation. 

In the case of body of revolution we can find the solutions of equation of surface 1)1.6)1 in the 
explicit form 

n = fi (73)71, r 2 = fi (73)72, r 3 = /2(73), (2.1) 

where ^(73), i = 1,2 are function subjected to the differential equation that defines the meridional 
section 

dfo = f l-7 3 2 ffi (22) 
d l3 11 73 d 73 ' 1 • ' 

If we denote the main central tensor of inertia as I = diag(ii, I\, 1$), (I\ = I2), then we can 
explicitly calculate the density of invariant measure of equations 1)1.1)1 . 1)1.21) . It exists for arbitrary 
functions 71(73), 72(73) that define the surface 

1 1 ^2 ^\ 

~ y/hh + m(r, Ir) ~ y/hh + hmf({l - 7a 2 ) + hmfi ' 

Remark. For the equations motion in variables u>,j the density of invariant measure differs from 1)2.3)1 
by factor det Iq, where Iq = I + ?n(r 2 E — r <g> r) is the tensor of inertia with respect to the point of contact. 
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In the case I\ = I2 we have detlQ = (Ji + mr 2 )(lil^ + m(r,Ir)). It is interesting that measure (|2.3[) contain 
one of these factors. 

It is easy to verify that under the above conditions the equations of motion also have the symmetry 
field v defined by the differential operator 

9=h - h A- hh A +ii A-'<^ (2A) 

It corresponds to the invariance of the system with respect to rotations about the axis of dynamical 
symmetry. Using this field we can reduce the order of system. For that we should choose the integrals 
of vector fields (|2.4|) as reduced variables to present the equations in the simplest form. After a number 
of tries, we choose the following reduced variables 

(M, r) f 2 
K x = V ; ' = M l7 i + M 272 + ^M 3 
/1 Ji 

K 2 = f = p(mhf 2 {Mvn + M 272 ) + (h + m/f )M 3 ) (2.5) 



M 27l - M i72 



V(l- 7 2 )(/ 1+ mr 2 ) 

In these variables the equations of motion of the reduced system have the following form 

73 = kK 3 
k 1 = -kK 3 plJl-(^f)')K 2 , 



k 2 = -kKzpmf 1 (h-f 2 )K 1 , 
^3 = - j2/2(1 fc _ 2)2 (/ 2 (/!(!- + 73/2) {mflKl + hK 2 2 )+ (2.6) 

+ l3 ffhKf - h{h{\ - 7 |) + 2 73 / 2 )^^+ 

+mffpf 2 (l - 7 |)( 73 /i/i - hh^K^ 



<9 7 

1-73 



where k = ,. 

y Ji + mr 

It is easy to show that these equations have the invariant measure with density p = k" 1 and 
integral of energy 



H = \(M,u) + U( l3 ) 

1 /V2 ^3 ^2 , ™/l /V ^2 \ 2 \ , 1^2 



(2.7) 



2/ 1 (l- 7 2 )V 1 m /2 2 ^ /x V 1 pm/1/2 

Moreover, for system (|2.6|) we have the following theorem 

Theorem. After the change of time kdt = dr vector field (|2.6|) become Hamiltonian 

^ = {xi,H}, x = ( 73 , Ki,K 2 ,K 3 ) 
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with Hamiltonian (j2.7j) and degenerate Poisson bracket: 

{73,^3} = 1, {^1,^3} : 

{K 2 ,K 3 } = -mpf 1 (f 1 -fl l )K 1 
(all the other brackets are equal to zero). 



J 3P(l-(f 



(21 



The proof of theorem is obtained by the direct verification of equations and of Jacobi identity. 
It turns out that equation ()2.6j) can be written in antisymmetric, almost Hamiltonian form (that 
sometimes is referred to as antigradient form) 



dxi 



J. 



dH 



'A > 



(2.9) 



where 














A 




-1 I 3 p[l- (^) )K 2 + Xu mph(h- &)K x + \v 



-hp[l-[fj )K 2 -Xu 
-mphtfi-ftiKx-Xv 




/ (2.10) 



11 



f?il(i-ii)K 3 



■K 2 ^--Ki 



m 



^-(h+mf^Kt ), 



h 

I?(l-T$)K 3 \fiP 

and A is an arbitrary function of (Ki, K 2 , K 3 ,^). At A = we obtain the degenerate tensor Jo 
corresponding to bracket (|2.8|) : although for A 7^ tensor (|2.1U|) is nondegenerate, it does not satisfy 
the Jacobi identity, i. e. it does not define a Poisson bracket. 
If A is chosen in the form 

A = a( 73 )^3, (2.11) 
then the tensor J = (A~ x Ja) satisfies the Jacobi identity, and the corresponding vector field 



v = (A" 1 J A )V J fi" 



(2.12) 



is Hamiltonian; at the same time the divergence of field (|2.12[) is nonzero. Thus, the considered 
nonholonomic system generate an example of Hamiltonian vector field with nontrivial measure p = 
= 0(73)^3. Note also that the function p = 0(73)^3 is a reducing multiplier in Chaplygin's terminol- 
ogy, and in this case it differs from invariant measure (|2.3|) . The close example of Poisson structure 
for the problem of ball's rolling on body of revolution was presented by Hermans |32j . 

Bracket (|2.8j) has two Casimir functions [H| that are integrals of motion; therefore, system (|2.6|) 
is integrable. The integrability and existence of linear integrals can be established by the different 
classical method: we divide the second and the third equation of system ()2.H|) by 73 and obtain the 
system of two linear non-autonomous first-order equations 



dKi 

dji 



-Phil 



dK 2 

G?73 



-mpfi(fi - f 2 )Ki 



(2.13) 



In somewhat different variables connected with semifixed axes, equations (|2.13|) were obtained by 
S. A. Chaplygin. Equations (|2.13|) do not contain the potential, which is presented only in the ex- 
pression of energy integral (|2.7|) . Using this equation we determine the dependence of nutation angle 
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on time after the solving of linear system (|2,13[) . In the general case this dependence has periodic 
oscillating character. 

Because the equations (|2.13j) are linear with respect to 73, the general solution can be obtained 
as the linear superposition 



(2) 



1,2, 



(2.14) 



where 5^(73), fff (73) are elementary solutions of (|2.13|) . Inverting expressions (|2,14j) with respect 
to Q we obtain the expressions for the lacking first integrals, which are expressed in the general case in 
terms of real analytic, but nonalgebraic (for example, hypergeometric) functions. Nevertheless, they 
are always linear with respect to Mj (i.e. with respect to generalized velocities). 

These integrals in some sense generalize the area integral (corresponding to the cyclic angle of 
precession) and the cyclic Lagrange integral (corresponding to the cyclic angle of proper rotation) [5, 
that exist in the classical problem of heavy symmetric top's motion about a fixed point (the Lagrange 
case). The presence of such integrals causes the great similarity of qualitative researches of these 
problems. 

Let's consider all known situations when integrals are algebraic or can be expressed through 
members of some known classes of special functions. 

Round disk (S. A. Chaplygin, P. Appell, D. Korteweg). Generally speaking, we consider 
a disk with the center of mass displaced along the axis of dynamical symmetry (Fig. |2J) • In this case 
functions (|2.1|) are explicitly expressed as 



h 



R 



h 



(2.15) 



where R is a radius of coin, a is the displacement of the center of mass along the axis of dynamical 
symmetry (Fig. EJ). 

The interesting fact in this case is the independence of mea- 
sure (|2.3|) from the phase variables p = const. For variables (|2.5|) we 
obtain the equations 



dK\ _ pmR' 



sm( 



K 2 , ^ = hp(sme + ±cose)K 1 , (2.16) 




where p = + I\mR? + I^ma 2 ) These two linear equations are 
reduced to one linear second-order equation with respect to 1V3 



Fig. 2 



- Ctg 0^3 + p 2 mR{R + a ctg 0)huj3 = 0. 

d6 z at) 



(2.17) 



At a = with the help of substitution cos 9 = 1 — 2x equation 1)2.17(1 is transformed to the hypergeo- 
metric type equation 

x(l - x)^ + (1 - 2x)^ - p%mR 2 u; 3 = 0. 
dx 2 dx 



In the papers |17l I24j the following result was shown: the disk for almost all initial conditions do 
not fall on the plane. In [2j the similar result was obtained for the nonintegrable problem of heavy 
disk's rolling on a slopping plane. 

The results concerning the stability of stationary motions and the qualitative analysis of motion 
see in the papers ^3 |H] and also in the book [IR- 
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Dynamically symmetric ball with the displaced center of mass (E. Routh, S. A. Chap- 
lygin). In this case 

h = R, f 2 = R 73 + a , (2.18) 

where R is the radius of ball, a is the distance from the center of mass to the geometrical center. The 
measure p is not constant any more 

p = (hh + hmR\l - 7 |) + I 3 m(R l3 + a) 2 )" 1/2 , (2.19) 

and the equations for K±, K 2 become trivial: K\ = 0, K 2 = 0, i. e. the expressions 

Ki = uzp^ 1 = p~ l (mR 2 j 3 (M, 7) + J1M3 + maR((M, 7) + M 373 ) + ma 2 M 3 ) = const, 

1 n (2-20) 

K 2 = -k(M, r) = (M, 7) + 75 M 3 = const. 

are integrals of motion. 

The integral K 2 represents the Jellett integral. This integral is also present under the arbitrary 
law of friction at the point of contact The integral K\ was found by E. Routh in 1884 |23| and its 
form is a little bit mysterious. It was also indicated by S. A. Chaplygin in the paper [2^]. Once again 
we shall note that both integrals are linear with respect to the velocities. They are the immediate 
generalizations of the cyclic integrals corresponding to the precession ip and to the proper rotation ip, 
but have no such natural dynamical origin. The integral K 2 sometimes is referred to as the Chaplygin 
integral. 

Remark. For axisymmetric bodies we can also indicate the other cases of existence of simple quadratic 
integral of the form 

F = aKf + bK%, a,b= const. 
Obviously, we have to require in addition to condition l|2.2|) the following one 



K^K 2 n*/i(/i-/a) 



= A = const. (2.21) 



The general solution of equations l|2.21|) and (|2.2fl is expressed in hypergeometric functions. Among the ax- 
isymmetric figures of the second order only the ball with displaced center satisfies these equations. The closed 
bounded curves satisfying (|2.21() and (|2.2() and different from ball arc similar to the ovals. Below we show that 
in this case the simple quadratic integral exists in the totaly symmetric case. (Ji ^ I2 7^ I3 h)- 

Three-dimensional point maps in nonholonomic mechanics. Before we consider the 
following cases of the body's motion, we shall present some general construction that let us to establish 
relations between equations (jl.2[) . H1.3|) to some point one-to-one map in three-dimensional space. 

We present the computer analysis of this map using the numerical integration of the indicated system 
at the fixed value of energy. Using this method we can find out and give a visual interpretation to 
various possibilities of existence of measure and integrals in their various combinations. 

To construct the three-dimensional map we use the Andoyer-Deprit variables (L,G,H,l,g,h), 
which were regularly used in our book for computer (and analytical) research of Euler-Poisson, 
Kirchhoff and other Hamiltonian equations. As against to nonholonomic situation, in the classical case 
these variables are canonical, and by virtue of the fact that the area integral is always present in the 
Euler-Poisson type equations, we can limit ourselves in the classical case to two-dimensional maps. 
The problems described above require two additional integrals of motion; therefore, it is necessary 
to use three-dimensional maps, and such maps are not necessarily possess an invariant measure (as 
against to Hamiltonian mechanics). Using the known formulas we make the transition from the 
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variables (M,-f) to the Andoyer-Deprit variables 

Mi = y/G 1 -L 2 sinZ, M 2 = \/G 2 - L 2 



cos / 



71 



72 



73 



§<»- 



§) + W 1_ (§) 2cos9 ) si ' 1 ' + \/ 1 -(f) 2si,19C08! ' 



G 



,g) (g 



a) + a*' 



cos g cos I — \ 1 



sin g sin/, 



1 



1 



cosg, 



(2.22) 



in which we can express the energy E = E(L, G, H, I, g). 

In the Euler-Poisson equations the value of H = (M, 7) is constant, but for equations (jl.lj) . I)1.2|) . 
1)1. 3jl this is not the case any more. We fix the level of energy E = Eq, then choose the intersecting 
plane, for example, as g = go = const, and obtain the three-dimensional map induced by sequential 
intersections of the phase trajectory with the chosen intersecting plane. We present the map in the 



^ 1, 



1. 



variables (L/G,H/G,l), because of its compactness by virtue of the fact that 

Typical examples of such three-dimensional maps are presented in Figs. [3 El 03 d It is obvious that 
because of the presence of one additional integral, the trajectories are situated 
on two-dimensional invariant manifolds of the point map, and the 
presence of two additional integrals imply the stratification of the 
three-dimensional space on invariant curves (Fig. |3J). In general 
case when both the integrals and the measure are absent, the com- 
plicated behavior of trajectories is possible. In this case random 
motions alternate with the asymptotic attracting properties typ- 
ical for dissipative systems. Note also that as against to (M,-f) 
the variables L,G,H,l,g,h are more convenient for the analysis 
of three-dimensional map, because in this variables the linear and 
angular components are separated, and they have the obvious ge- 
ometrical meaning (see 

One of examples of three-dimensional map is the well-known 
Smale- Williamson map. It does not preserve the measure, but is 
expressed by analytical formulas. Other examples can be obtained 
by the study of general (nonconservative) perturbations of two- 
degree Hamiltonian systems. 




Fig. 3. 



Rolling of balanced, dynamically nonsymmetric ball 
(Chaplygin ball [25J ) . The equations of motion of dynamically 
nonsymmetric ball with the center of mass coinciding with the 
geometrical center can be written in the form 



The three-dimensional map 
described in subsection [21 for the case 
of Chaplygin ball. The figure shows 
very clearly that all trajectories are 
situated on joint level surfaces of two 
integrals H = const and M 2 = const 



(h 



1, h = 2, h 



3) 



M 



dU 

M x u + 7 x — — , 7 = 7x0;, 
cry 

M = lu + D-y x (u x 7), D 



(2.23) 



ma 



where I = diag(ii, I2, 13) is the central tensor of inertia, U = U{"f) is the potential energy. Equa- 
tions 1)2.23)1 always have the measure with density p and the first integrals of the form 



P 



1 



VI -D(y, A 7 )' 
1 



(I + DE) 



E 



J i3 II 1 



(2.24) 



H = ±(M,u) + U(j), Fi= 7 ^ = l F 2 = (M, 7). 
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At U = there exists the additional integral F% = M 2 and the problem become integrable 
(S. A. Chaplygin, 1903 US]); the corresponding three-dimensional map is presented in Fig. 01 
It was shown in the paper that this problem is still integrable with a Brun potential 



U = ife(I 7 , 7)- 



The integral F3 in this case has the form 



The authors indicated in [5] that for any potential U the change of time dr = pdt in equa- 
tions (|2.23|) makes them Hamiltonian with a Poisson bracket, which is nonlinear with respect to the 
phase variables {M,~f) and has the form 

{Mi,Mj} = e ijk p~ 1 (M k - gj k ), {Mi, 7j} = £ijkP~ l lk, {li, 7j} = 0, 

D( 7 ,AM) (2.25) 



g = D(w,7) 



1-D(7,A 7 )" 



Bracket (j2.25|) is degenerated; its Casimir functions are integrals F\,F2 (|2.24|) . The Hamiltonian 
corresponding to bracket (|2.25|) is obtained from energy l|2.24[) expressed as a function of the moments 
by the formula 

H = \{M,AM) + \g(AM, 1 ) + U( 1 ). (2.26) 

After the change of variables K = pM the Poisson bracket and the Hamiltonian are represented in 
the form 

{Ki,Kj} = e ijk {K k - Dp 2 (K,~f)a k j k ), {K i ,j j } = e ijk j k , { r y%,Jj} = 0, 
H = \p' 2 {K,AK) + \D(XK , 7 ) 2 + U(j). 

Thus, at the zero level (if, 7) = bracket (|2.27f) passes to the bracket described by algebra e(3). on 
which we can write the Euler - Poisson and Kirchhoff equations [H] ■ 

Note that for the considered problem the density of measure p is the reducing multiplier (by 
Chaplygin [221). With its help the nonholonomic equations are reduced to the Hamiltonian system. 
Chaplygin himself used such reduction integrating the equations of motion of nonsymmetric ball; 
as a preliminary he introduced a nonholonomic analog of spherocon variables. It is possible to do 
these operations in inverse order [221: one makes at first the change of time dr = pdt to receive a 
Hamiltonian system, and then introduces the usual spherocon variables and use the Hamilton- Jacobi 
method. 

As against to Poisson structure Q2.8JI related to the system reduced on the field of a symmetry 
corresponding to proper rotation, structure (|2.12j) is related to complete system (jl.ip . (|1.2f) . Un- 
fortunately, we were unable to generalize (to lift) reduced structure Q2.8JI to such complete system. 
Possibly it is either too difficult or some dynamic effects prevent such generalization. Unfortunately, 
the dynamic effects preventing the reduction to Hamiltonian form are very purely investigated [1]. 

Rolling of unbalanced, dynamically nonsymmetric ball on plane. In this case equations 
(jl.2() can be written in the following convenient form 

{M = M x u + mr x (u x r), 
r = r x u — axw = (r- a) x u (2.28) 
M = lu + mr x (u> x r), 
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where a is the vector connecting the center of mass with the geometrical 
center r = Rry + a (see Fig. |1J). It turns out that in the case a / the 
integral F$ = M 2 of system (j2.23j) can be directly generalized and written 
in the form 



mr 2 (M 1 uj) = M 2 - 2mr 2 H, 



(2.29) 



where H 



-(Af,u;) is the energy integral. Though this integral is simple 




TTTTTTTTmTrTTTTTTTTTTTTT 



enough, but probably it was not known earlier. Note only that under the Fig. 4 

additional requirement of the dynamical symmetry the Jellett integral and the Chaplygin integrals 
were found (see subsection |2J). Integral (|2.29f) can be considered as their generalization for dynamically 
nonsymmetric situation. 




Fig. 5. One of trajectories in the problem of rolling of unbalanced ball on plane. The figure shows clearly that 
all points are situated on some surface; the condensations of points correspond to asymptotic approximations 
of the trajectory to periodic solutions. The trajectory goes out from the top and approaches to the three points 
in lower part of surface 




Fig. 6. Three trajectories in the problem of rolling of unbalanced ball on plane. The figure shows very clearly that 
points are situated on two-dimensional surfaces (corresponding to the level of integral |2.2Hjl ). The condensation 
of points corresponds to an asymptotic approximation to some periodic solution. 

We were unable to obtain generalizations of this integral for cases with gyrostat and Brun field. 
Note also that for a ^ there is probably no measure. This is illustrated in Figs. |21 El These 
figures shows the asymptotic trajectories of the point map situated on the two-dimensional surface of 
integral (LT291 . 

An arbitrary body with a spherical central ellipsoid of inertia. I\ = I2 = I3 = /U, 

fi = const. 
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This problem in general case requires both additional integrals, nevertheless, there is always an 
invariant measure. For the first time this fact was noted by V. A. Yaroschuk [2S]. In this case the 
equations written in variables u>,-y are more convenient. They have form 

(u + mr 2 )w = m(r + u x r)(r, u) — muj(r, r) + ■y x — — , 

5 7 (2.30) 

7 = 7x01. 

Equations (j2.3Uf) have an invariant relation (cj, r) = 0, which is used for simplification of some 
calculations. 

The density of invariant measure for these equations indicated in |2H] can be represented in the 

form 

p={p + r 2 )^. (2.31) 



We were unable to reproduce the generalization of this measure for the case of body's rolling on sphere 
indicated in the paper [2Hj. Probably, this generalization does not exist. Also it is not known, what 
nontrivial integrable cases can be obtained with the help of measure (|2.31|) . and whether the system 
has any Hamiltonian origin, possibly after the appropriate change of time. 



Gyrostatic generalizations. Following mainly the paper by S. A. Chaplygin |26j . we present 
the generalizations of the indicated problems for the case with additional uniformly rotating balanced 
rotor. The corresponding system can be interpreted as a nonholonomic gyrostat. The gyrostatic 
effect can be also obtained by an addition of multiply connected cavities completely filled by the ideal 
incompressible liquid possessing nonzero circulation into the body 0. In the described case, equation 
for the moment (|l.lj) can be presented as 

M = (M + S)xw + mfx(uxr) + M Q , 

where S is the constant three-dimensional vector of gyrostatic moment. It is easy to verify that the 
addition of rotor does not influence the existence of invariant measure with the density depending on 
the positional variables 7. 

a) Body of revolution. The equations of type (|2.13|) for the rotor with gyrostatic moment S = 
= (0, 0, s) directed along the axis of revolution in variables (|2.5j) have the form 

g = -/3.(l-(f)> 2 - S , d -^ = -m Ph ({ h -f 2 )K l + h s). (2.32) 

Equations (|2.32|) were obtained in less convenient form by S. A. Chaplygin [2^] . The density of invariant 
measure is also defined by equation ()2.3|) . 

Let's consider sequentially the generalizations of the indicated earlier cases of disk, ellipsoid, and 
ball with the displaced center. 

b) Round disk. Now equation (|2.17|) have to the following form 

^ - ctg + mRh{R + a ctg 6) P 2 lj 3 = smRp 2 (R + a ctg 9), 
d9 z at) (2.33) 

p = (hh + hmR 2 + I 3 ma 2 )~^, 

and at a = in general case it is reduced to non-homogeneous (for hyper geometric equation. 

c) Ball with the displaced center of mass. Here system ()2.32(l has the form 

^ = -s, ^ = -mpR(R 73 + a)s, (2.34) 
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where p is defined by relation (|2.19|) . We can immediately show the first integral that generalizes the 
Jellett integral 

F = K\ + 573 = const. (2.35) 
The second integral generalizing Routh (Chaplygin) integral has more complicated nonalgebraic form 



(h - I^p- 1 ^ - s jp" 1 - h^j arctg ^ 1 -JIl 1 -p(Rr Y3 (I 1 - J 3 ) - al 3 )j j = const. (2.36) 

In integral (|2.36|) we assume 1\ > I3. For I\ < I3 the integral contain hyperbolic functions. Inte- 
gral (|2.36f) was explicitly presented by A. S.Kuleshov [T^]. It is essentially simplified at I\ = 1% = p 
and has the form 

p- 1 (3/ia* - S + nR 2 -*™ 2 -™^ ) = const . (2 . 37) 
V ma ' 

The form of integral is even simpler for the case of balanced homogeneous ball (a = 0): 

UJ3 + 7j — ^^^=^^^ = const . 
^ \J p(p + mR 2 ) 

This simple integrable generalization was indicated by D. K. Bobylev jH] (some additional simpli- 
fications in the case of explicit integration were also indicated by N. E. Zhukovsky |12j). 

d) Dynamically nonsymmetric ball. The most general gyrostatic generalization for the case of 
Chaplygin ball was suggested by A. P. Markeev [^Oj. The equations of motion and the integrals are 

M = (M + S) x u, 7 = 7x0;, 

1 , (2-38) 

fl-=±(Af», F 1=7 2 = l, F 2 = (M + S,M + S), F 3 = (M + S, 7 ), 

where S is the constant three-dimensional vector of gyrostatic moment. Note that we were unable to 
generalize Poisson structure (|2.25|) to system (|2.38|) for S / 0. 

Rolling of ellipsoid on plane. It turns out that in the problem of rolling of balanced ellipsoid, 
which axes are principal axes of inertia, there are cases of existence of specific invariant measure, that 
are defined by restrictions on ratios of moments of inertia and semiaxes of the ellipsoid of the surface. 
This measure has found by V. A. Yaroschuk in j^Hj ■ 

For the problem of ellipsoid's rolling, two cases of existence of the invariant measure were already 
indicated in subsections |2 and [3 They are accordingly measures of balanced, dynamically nonsym- 
metric ball ()2.24|) and of arbitrary body with the spherical central ellipsoid of inertia ()2.31|) . It is 
interesting that no obstacles to existence of the analytical invariant measure of the general rolling 
problem of balanced ellipsoid, which principal axes are principal axes of inertia, are not found yet (as 
against to the case of Celtic stone JS1). It is possible, that this measure exists (at least in this situation 
there is no asymptotical behavior typical for the Celtic stones), but is complicated and nonalgebraic. 

For the surface of ellipsoid (r, B _1 r) = 1, where B = diag(&i, 62 ; ^3)) h are squares of larger 
semiaxes, we have the explicit expression 

Bj (2.39) 



V( B 7, 7) 

If the cental tensor of inertia has the form 

I = pE + AmB, (2.40) 
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Fig. 7. Some trajectories in the problem of rolling of ellipsoid with the spherical tensor of inertia on plane. The 
figure shows clearly the regular trajectories filling some curves that enclose stable periodic solutions (permanent 
rotations), which in this case are degenerated as it was stated above. A random layer (which is obtained from 
one trajectory) in this case is not situated on any surface. 

then equations ()1.2|) have a measure only in the case A = and A = 1 with the density (in 

variables M, 7) 



Remark. For the case A = 0, measure H2.41|l was already indicated in subsection It is defined by 
formula (j2.31|l and is present for any surface of the body. The differences in powers of expressions l|2.41l) 
and H2.31fl are connected with the various systems of variables (M,~y) and (^,7) and with the corresponding 
transformations of densities of invariant measures. 

Note that if equality ()2.4(Jj) is fulfilled, then the motion equations have a two-parameter set of verti- 
cal permanent rotations at arbitrary real A (in other cases this set is one-parameter). A. V. Karapetyan 
in showed that the conditions of existence of such sets are even a little bit wider and have the 
form 



where I = diag(/i, I2, In addition to tensor (|2.4U|) conditions ()2.42|) are fulfilled in the case of non- 
holonomic Chaplygin ball (61 = 62 = ^3 = R 2 ) when there exist both the measure and integral ()2.24j) . 
and also the conditions hold in the axially symmetric situation I\ = I2, b\ = 62- Unexpectedly, 
equality (|2.4(H) is also the necessary (but, generally speaking, insufficient) condition of the integrabil- 
ity of equations of motion for the case of ellipsoid on ideally smooth plane [S]. Note the interesting 
fact that for A / and A / 1, both the measure and obstacles to its existence are not found. The 
three-dimensional section in the case A = is presented on Fig. [7| 

All the results described above that are connected to the rolling of body on plane are presented 
in table 1. 

3. Body on a Sphere 

Now we consider systematically the situations analogous to one in the case of plane that originate in 
the problem of rolling of body on a sphere. First of all we shall note that kinematic equation 
can be written as 




(2.41) 




(2.42) 



ijk 



7 = 7 x u T kr 



k = 1/Ro 



(3-1) 
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Table 1. Rolling of body on plane 



tensor of inertia 


dynamically nonsymmetric case 
h^hlt h + h 


axial dynamical symmetry 
h =h,U = U(j 3 ) 


total dynamical 

symmetry 
h= h= h = tt 


surface 
of body 


ball 


ellipsoid 


an arbitrary body of 
revolution 


round disk with 
sharp edge 


unbalanced ball 


arbitrary 


geometrical 
and dynamical 
restrictions 


the center of mass 
coincides with the 
geometrical center 


the center of 
mass does not 

coincide with 
the geometrical 
center 


the axes ot 
dynamical and 
geometrical ellipsoid 
coincide 
I = + mB 


the geometrical and dynamical axes coincide 
and contain the center of mass 




measure 


(l- J D( 7 ,A 7 ))- 1 /2 


unknown 


(fi + mr 2 )~ k 


(hh + m(r,Ir))- 1/ 2 


const 


(Iih + mirjr))- 1 ^ 


(/j, + mr 2 )'^ (for 
variables u> and 7) 


additional inte- 
grals 


M 2 = const 

(M,7) = const 
(two integrals) 


M 2 - mr 2 (M,Lj) = 

= const 
(one integral) 


none of integrals are 
found 


two integrals are 
obtained from the 
solution of system of 
two linear 
equations 12.131 


two integrals are 
obtained from the 
solution of 
hypergeometric 
equation 12.1711 


uJa/p = const 
(M, r) = const 


none of integrals are 
found 


integrable 
addition 
of gyrostat 


possible 
(A. P. Markeev 
1986) 


it seems to be 
impossible 


the measure is 
preserved 


S. A. Chaplygin 
(1897) 


S. A. Chaplygin 
(1897) 


at 1\ — Li = 73 the 
gyrostat was added 
by D.K. Bobylev, 
at I 1 = I 2 =£ I 3 by 
A. S. Kuleshov 
(2000) 


the measure is 
preserved 


Hamiltonian 
form 


the system is 
Hamiltonian after the 
change of time 
(A. V. Borisov, 
I.S.Mamaev, 2001) 


it seems that the 
system is not 
Hamiltonian 


unknown 


the reduced system is Hamiltonian after the 
change of time, defined by the reducing 
multiplier 
(A. V. Borisov, I.S.Mamaev, 2001) 


unknown 


authors 


S. A. Chaplygin 
(1903) 


A. V. Borisov, 
I. S. Mamaev 
(2001) 


V. A. Yaroschuk 
(1995) 


S. A. Chaplygin 
(1897) 


S. A. Chaplygin 
(1897), P. Appell, 
D.Korteweg (1898) 


E. J. Routh (1884), 
S. A. Chaplygin 
(1897) 


V. A. Yaroschuk 
(1992.) 


generalizations 
and remarks 


the integrable addition of 
Brun field is possible 
(V.V.Kozlov, 1985) 
the Hamiltonian form is 
preserved for arbitrary 
fields with the loss of one 
integral 


the Brun field 
can not be added 
(preserving the 
integral) 













H 
H 

PS 

O 

fi 
fi 

o 
o 

H 
O 

O 

> 

s 

o 
a 

a 
o 
a 

* 

o 

> 

> 

> 

a 

> 

*d 
K 

5 



PS 
> 
Pi 

o 

o 
a 

> 



Remark. The cases of existence of the corresponding (tensor) invariants are indicated by gray color in the table. The partial filling corresponds to the uncomplete 
set of integrals. 



A. V. BORISOV, I. S. MAMAEV 



where the "minus" sign means the rolling on interior of surface of sphere (Fig. IHla), and the "plus 
sign means the rolling on exterior of surface of sphere (Fig. |H]b). 




Fig. 8. 



Rolling of body of revolution. The body of revolution is defined analogously by formu- 
las (|'2.1[) , and we can obtain the explicit form of the density of invariant measure 



Po 



K73) = Po(l - khf(l - kf 2 ), 
[hh + m(r,Ir))~ 1/2 = (hl 3 + mlj^l - 7 I) + ml 3 fhiy 1/2 



(3.2) 



and establish the presence of the symmetry field v with operator ([2.4)1 . 
Reduced system (j2.13j) in the variables of type (|2,5|) 



K x = (M,r)(l-fc/ 1 )/r 1 , K 2 



U3P0 



(3.3) 



has the form 



1 - kfi dK 2 
Po ^73 



1 dKi 

Po d7 3 



hi 1 



+ kf 1 {h-h)(l-kf 2 ) 



-m 



- f 2 ~ kfiffiKx + mkpohf! 73/2 + V 1 " ^ /1V 1 " ^ 



(3.4) 



Equations ()3.4|) in somewhat different variables connected with the semifixed axes were obtained by 
P. V. Woronetz who generalized the Chaplygin arguments to the case of sphere. 

Integrating the equations of system (|3.4[) we determine the dependence on time of the nutation 
angle 9 = arccos(73) using the quadrature of the energy integral 

H = i^+^7l+ 



2 1 



+ - 



2h{l- 1 i)\(l-kf l Y mf( 



Jo rat 
3 Kl+ J - 



7i 



K 2 



h vl-fe/i mp/1/2 



+ C/(73). 



For the problem of rolling of the body of revolution on sphere we can describe the reduced system 
in variables ()2.5|) and the corresponding Poisson structure similarly to the case of plane. We do not 
present the corresponding calculations because of their bulkiness. 

Similarly to the problem of rolling on a plane we shall consider some special cases. 
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The problem of rolling of round disk with the center of masses displaced along the axis 
of symmetry in general case is not reduced to the hypergeometric equation any more, but the system 
nevertheless becomes simpler; thus, the measure is not constant any more. Same as above, /i,/2 are 
defined by relations (|2.15|) . and the density of invariant measure has the following form 

p = p (l kR - ) , p = (l x I 3 + m(hR 2 + /3a 2 )) / = const, 

V vi-73/ v ' 

where R is the disk radius and a is the displacement of the center of mass. 
The second-order equation has the form 

+ ctg 9 1 + -7— - a uj 3 — [R + a ctg 9 + — w 3 = 0, (3.5) 



d# 2 V sm9-kRJ -. kR \ h sin I 



sm( 



where 73 = cos 6. 

This equation in the particular case I 3 = 27i, a = (the homogeneous balanced disk) was 
obtained by P. Woronetz [3U] , and at a = it was investigated in [33] by methods of the qualitative 
analysis. In particular, the stability of stationary motions and the probability of falling of disk on 
sphere was investigated and this probability is found to be equal to zero. In the paper [31] it was 
also shown that, as against to the nonholonomic problem, the the Hamilton system describing the 
motion of disk on a sphere with absolute (ideal) sliding that seems to be simpler is not integrable any 
more, and its behavior has random properties. The number of degrees of freedom for this system is 
increased in comparison with the case of plane, where the indicated system is integrable because of 
the preservation of impulse's horizontal component, and this is the reason for such result. 

Ball with displaced center. Here functions f±, f 2 are also defined by relations (|2.18|) . and the 

expression for measure p is the same as in the case of plane: 

p = (liJ 3 + hmR 2 (l - 7 |) + I 3 m(R 73 + a) 2 ) 
In variables ()3.3|) equations ()3.4|) have the form 

±K[ = -kR(h-I 3 )(l-kR)K 2 , 1^ = ^^^. (3.6) 

With the help of these equations we obtain two linear with respect to K\ , K2 nonalgebraic integrals 
of the form 

F 2 = ( y/ m (I 3 - h) K 2 + K x )e x \ F 3 =^m(I 3 -h)K 2 -- [ ^ n K i y~^, (3.7) 

where A 2 = mk 2 R 4 (I 3 — 1±), r = f ^0(73) ^73, and the additional quadratic algebraic integral (depen- 
dent on ^2,^3) 

F = F 2 F 3 = ™** K* + (h - I 3 )K 2 2 . (3.8) 
(1 - kR) z 

The integrals F 2 , F 3 are new and generalize Routh and Jellett integrals (|2.2U|) . Integral (|3.8|) was 
found by A. S. Kuleshov [18| . In the case of spherical tensor of inertia ^3 = 1\ we have 

K x = (M, r) = Mi 71 + M272 + M 3 ( 73 + -|) = const, 

, (3.9) 
(1 - kR)K 2 - kmR 3 K 1 / ,9(73) dj 3 = const, 



i. e. K\ coincide with the classical Jellett integral. 
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Balanced, dynamically nonsymmetric ball on sphere. 

We shortly describe an integrable case connected with the problem 
of rolling of balanced, dynamically nonsymmetric ball on sphere anal- 
ogous to the motion of Chaplygin ball on plane (see ()2.23j) . ()2.24(1 ). 
This system is defined by the equations 

M = M x cj, 7 = kj x uj, 
M = (I + D)u- D~r(u,~f), 

R 



(3.10) 



where k 



R 



a is the radius of ball, m is its mass, D = ma 2 , 



Fig. 9 



R is the radius of the fixed sphere (see Fig. (In the case of plane 
R — > oo and k = 1.) Equations (|3.1U|) always have the integrals 

1 



H 



[M, u>) 



1. 



M l 



and the invariant measure with density p ()2.24j) . We need one more integral for the complete integrabil- 
ity of system (such as the area integral at k = 1). It exists only under the additional condition a = 2R 
found by A. V. Borisov in [7j, which corresponds to the rolling of balanced, dynamically nonsymmetric 
ball on the interior of the fixed sphere (see Fig. |§J). This integral have the form 

F = (M,A 7 ), A = E- 2(tr(I + J D))- 1 (I + L>) (3.11) 

and can be generalized to case (|3.1U|) with the additional Brun potential U = ^(17,7) [7]. Using 

the transformations M = AM, 7 = 7 under condition a = 2R we transform equations (|3.1U|) 
into the equations describing the motion of Chaplygin ball on the horizontal plane. For arbitrary 
parameter k and potential U the Poisson structure of equations (|3.1U|) similar to one indicated at k = 
= 1 in subsection 20 (formula (|2.25^ is still unknown. 

Unbalanced, dynamically nonsymmetric ball on sphere. In this case there exists one 
(and only one) quadratic integral generalizing the corresponding result on plane. The equations are 

M = M x uj + mr x (u x r), 

1 ^. ( 3 - 12 ) 



where r = R~y + a, M 



Icj + mr x (u x r). It has exactly the same form 
F = M 2 



mr 2 (M,u). (3.13) 

Note also that for a 7^ it seems that there is no invariant measure, and also generalizations for 
the cases with additional gyrostat or with the field of Brun problem are unknown. In the case of 
rotational symmetry I\ = I2, integral (|3.13|) was indicated by A. S.Kuleshov and the acquaintance 
with this result has induced authors to the analysis of dynamically nonsymmetric situation. 

Gyrostatic generalizations. Let's briefly discuss the integrable gyrostatic generalizations. So 
to preserve the integrability for the body of revolution's case we should direct the balanced rotor with 
the moment S along the axis of dynamical symmetry. In variables (|3.3fl the analog of system ()2.32j) 
has the following form 

- ) I +kh(h-h)(l-kf 2 ) 



Po Kl 



kfi 



-m 



v 3)) +kh{h 

+sp- 1 (l-kf 2 ){l-kh), 
kmpohf? 73/2 + V l- 7 3 2 (/iyi^ 
- /a " kfxf^Kx - smfMl - kf' 2 ){l 



K' 



(3.14) 



K 2 - 



kfi 
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where po is defined by relation (|3.2j) : thus, the linear system is non- homogeneous 

In the case of round disk, the integrals and equations (|3.14[) can 
not be essentially simplified. In the case of ball with the displaced cen- 
ter, fi , f 2 are defined by equations H2.19j) , and equations (|3.L4f) are 



Po 1 



K' 



-kR(I x -h){l-kR)K 2 
mR 5 k 



J_ 

Po 



R(l - kR) 2 s, 



(3.15) 



1 - kR 



Ki - mR(l - kR){Rj 3 + a)s. 



The first integrals in this case can not be written with the help of ele- 
mentary functions. 

In the case of total dynamical symmetry I\ = 1% with the help of 
equations (|3.15j) we obtain the explicit integral of Jellett type 



F 2 = K l + {l-kRf 73+ 



R 



{MR^ + Sil-kR), r) (1-kR) = const. 




Expressing from this integral K\ and substituting the expression in the 
first equation of system (|3,15|) . we obtain the explicit quadrature for 
^1(73). The gyrostatic generalizations of integral (|3.11j) for the case of 
nonsymmetric ball are unknown. 

Rolling of body with partially flat surface on sphere. Let's con- 
sider one more problem connected with the rolling of body on a sphere, which 
has no analog in the case of plane. We consider the rolling of the body with 
the flat foundation on the exterior surface of sphere (Fig. Illj) . This problem for 
the first time was considered by P. Woronetz, who indicated the integr ability 
of the problem in the case of rotational symmetry. Let's write the equations 
of motion in the frame of references connected with the body in the case when 
the field of force is absent: 



Fig. 10 



M = M x u + mr x (a; X r), M = Iuj + my x (u> X 7), 



(3.16) 



where I is the central tensor of inertia, m is the mass of body. 

For the flat part of surface we have r = (n, r 2 , z), z = const, 7 = (0, 0, 1), 
where n, r 2 are the projections of position of the center of mass on the flat 
foundation (see Fig. 111(1 . and using the fact that 7 = we obtain for them 
using equation (jl.3|) 

ri = k~ 1 uj 2 , f 2 = —k~ l uj\, k = R~ . 




Fig. 11 



(3.17) 



We can specify two cases, when equations (|3.16|) . (|3.17(l have invariant measure. 

a) z = — the center of mass is situated on the contact plane, and tensor of inertia I 
= diag(ii, I 2 , I3) is arbitrary. For the variables M, r the density of invariant measure is 



— 1 /2 

p(M,n,r 2 ) = (iz + mr^j 
It can be also written for the equations in variables u, r%, r 2 : 

p{u,n,r 2 ) = {hl 2 + m(r,lr))(I 3 + mr 2 )^ . 



(3.18) 



(3.19) 



Remark. In this form under additional and not essential restriction I3 = Ii + I2 (i- c. in the case of flat 
plate) the invariant measure was indicated by V. A. Yaroschuk |28j . 
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In the considered case in addition to the measure, equations (|3. 16f) . (|3. 17|) have one (and only 
one for I\ 7^ I 2 ) first integral 

F = M 2 — 2mr 2 H, H = |(M, u), (3.20) 

independent of energy integral. It was not known earlier. 

Besides, it turns out that integral (13.20(1 is transferred without modifications to the case z 7^ 0, for 
which the measure is not known (and probably does not exist). Note that we already see the similar 
situation in the problem of rolling of unbalanced, dynamically nonsymmetric ball (see subsection E| . 
The lack of measure for z 7^ obviously prevent the existence of Hamiltonian form. For z / though 
the measure exists, the system is probably also not Hamiltonian even after the appropriate change of 
time. 

b) I\ = I2, z 7^ 0, i. e. the center of mass is situated on the axis of dynamical symmetry. In 
variables (M,r) the density of invariant measure is 

p= (I 1 I 3 + ro(r,Ir))-V (3.21) 
This system is already integrable. Indeed, in variables 

K x = (M, r) = Mm + M 2 r 2 + M 3 r 3 , 

to 3 mz{M 1 r 1 + M 2 r 2 ) + (h + mz 2 )M 3 

2 = ~7r = — 

p y/hh + m{r,lr) 

the equations of motion of system have linear form 

dK 1 = W3 R dK2 = m Ri {3 22) 

du 2R^J h(h + mz 2 ) + hmu du 2R^I 3 (h + mz 2 ) + hmu 

P. Woronetz in the paper |38| noted that their explicit solution can be obtained in elementary functions. 
This solution defines two linear with respect to M additional first integrals. 
There are two different methods of solving of (|3.22|) . 

1) ii < / 3 - 



K\ = \J 1\ — h(—c\ cos ip(u) + c 2 sin (f(u)) , K 2 = \/m(ci sin <p(u) + c 2 cos (f(u)) , 

2 (h - I 3 )(hl 3 + I 3 mz 2 + hmu) (3.23) 
<p= , Cl ,c 2 = const. 

2) I\ > ^3, In this case we should use hyperbolical instead of trigonometrical functions 
In this case there is also simple, but dependent quadratic integral 

F = mK 2 + {h - I 3 )K 2 , (3.24) 
that evidently is a particular case of (|3.20|) . The energy integral can be presented in the form 

H = i(M,w) + U(r) = |^ - 2 u 2 + 

K U (3.25) 

+Ar({h + mz 2 )K\ + (hu + I 3 z 2 )Kl - -%-Vhh + hmz 2 + hmuK x K 2 ) + U{u), 
2I(u V l{u ' 
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and using this form we obtain after the integration of system 1)3.22(1 the explicit quadrature for u. 

In the axisymmetric case we can add rotor with the moment S = (0, 0, s) along the axis of 
dynamical symmetry; measure (|3.21() is preserved, and equations (|3.22[) are 

2R^ = k^M K 2 -s, 2R^ = rn(K 1 + zs) 

du y/I-sih+mz 2 ) + hmu du h{h + mz 2 ) + hmu 

The general solution of system (|3.2ti|) is presented as a superposition of the solution of homogeneous 
equation at s = (|3.23j) and the particular solution of non-homogeneous equation (|3.26|) of the 
following form 

(part) _ Sz(h + Ji(l + R/z)) (part) _ S ^ h(h + mz 2 ) + hmu 

Thus quadratic integral (|3.24|) in this case is presented as 

F = m{K x - K[ part) ) 2 + (h - h){K 2 - K { 2 part) f. 

Remark. We have noted in the paper two nontrivial cases of non-existence of one quadratic integral. The 
similar integral is present in the case of homogeneous ball's motion on the triaxial ellipsoid. The origin of such 
integrals is connected with the presence of similar integrals in the axisymmetric situation, when the system is 
completely intcgrable, and there are two linear integrals. However, in the general case the dependence on the 
positional variables in these integrals is complicated and is expressed in special functions. In three cases that 
we have found, the linear integrals are expressed in elementary functions and generate the quadratic integral, 
which is represented by rational function. This integral can be also generalized to the nonsymmetric situation. 



4. Conclusions 

In this paper we collect all the known at present cases of existence of invariant measure, integrals, 
Poisson structure for the equations of nonholonomic rolling of rigid body on plane and sphere. In all 
conceivable assortment of situations we have not found out any case when there are two integrals, but 
no measure. Possibly it is connected with the specificity of the equations of nonholonomic mechanics. 

Depending on the presence of this or that set of invariants there are qualitative distinctions in the 
behavior of system. The system can exhibit both typically Hamiltonian properties and the properties 
of conservative systems realized on the example of Celtic stones' problem. For the analysis of the 
indicated situations we use the method of three-dimensional Poincare maps and exactly this method 
at first let us to find out integrals (|2.29|) . I|3.13|) . (13.20(1 numerically, and then to obtain their explicit 
form. Undoubtedly, the research of three-dimensional point maps in the cases of presence and lack of 
the measure both from analytical and from the computational point of view allows to find out many 
remarkable effects in nonholonomic systems. 

From the point of view of this approach the problem of global evolution of the Celtic stone jl4H 21j 
is the most interesting. In this problem we have to study the invariant and asymptotic manifolds on 
the level surface of three-dimensional maps (not preserving an area) realizing under various restrictions 
on the parameters of system. 

Authors thank A. V. Karapetyan for useful discussions and A. A. Kilin for the help with realization 
of numerical experiments. 
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Table 2. Rolling of the body on sphere 



tensor of 
inertia 


dynamically nonsymmetric case Ji 7^ I2 7^ h 7^ Ii 


dynamical symmetry X\ = I2, U — £7(73) 


surface 

of the body 


ball 


body with partially flat surface 


an arbitrary body of 
revolution 


round disk 


unbalanced ball 


body with partially 
flat surface 


geometrical 
and dynamical 
restrictions 


the center of mass coincides with 
the geometrical center 


the center of 
mass does not 
coincide with 

the 
geometrical 
center 


the center of 

mass is 
situated on 
the contact 

plane [z — 0) 


the center 
of mass is 

not 
situated 
on the 
contact 
plane 


geometrical and dynamical axes coincide and contain the center of 
mass 




a = 2R 


measure 


(1-_D( 7 , A 7 ))" 1/2 
(V.A. Yaroschuk, 1992) 


unknown 


(V. A. Yaroschuk, 
1992, 
A. V. Borisov, 
I. S. Mamaev, 
2001) 


unknown 


(i-fc/i) 3 (i-fc0sA)') 

y/hh + m(r,Ir) 
(A. V. Borisov, I.S. Mamaev, 2001) 


(/i/ 3 + m(r,Ir))- 1/2 


integrals 


M 2 = const 
one integral 


M — const, 
(AT, A7) - 
— const 
(two integrals) 
(A. V. Borisov, 
1994) 


M 2 -mr 2 (M,w) - 
— const 
one integral 
(A. V. Borisov, 
I. S. Mamaev, 2001) 


M 2 - mr 2 (M,Lj) = const 
(A. V. Borisov, I.S. Mamaev, 2001) 


two integrals are 
obtained from the 
solution of system of 
two linear 
equations p.4t 
(P. V. Woronetz, 
1909) 


two integrals are 
obtained from the 
solution of second 
order equation $3.b\ 
(P. V. Woronetz, 
1909) 


two integrals are 
expressed in 
elementary 
functions l|M.7fl 
(A.S.Kuleshov, 
2000, A. V. Borisov, 
I.S. Mamaev, 2001) 


two integrals are 
expressed in 
elementary 
functions |J3.23^ 
(P. V. Woronetz, 
1911) 


integrable 
addition 
of gyrostat 


possible 
without loss of 
integral and 
measure 


not 
found 


not 
found 


not found 


not found 


possible along the axis of dynamical symmetry 


Hamiltonian form 


probably, the 
system is not 
Hamiltonian 


the system is 
Hamiltonian 
after the 
change of 
time 
(A. V. Borisov, 
I. S. Mamaev, 
2000) 


piUUaUlV. Lilt; 

system is not 
Hamiltonian 


probably, the system is not Hamiltonian 


the reduced system 


is Hamiltonian after 1 
multiplier (A.V.Boris 


he change of time defined by the reducing 
3V, I.S. Mamaev, 2001) 


generalizations 
and remarks 


the addition of 
Brun field is 
possible 
preserving one 
integral and 
measure 


the integrable 
addition of 
Brun field is 

possible 
(A. V. Borisov, 
Yu. N. Fedorov 
(1994)) 


the Brun field 
could not be 
added 
(preserving the 
integral) 


V.A. Yaroschuk 

found the 
measure under 
additional 
nonessential 
restriction 
h =h + h 








A. S. Kuleshov 
found only one 
quadratic 
integral 
dependent on 
two present 
linear integrals 


P. V . Woronetz 
found the solution 
in quadratures 
under additional 
nonessential 
restriction 

= 1/3, 

z = n 



> 
'< 

td 
O 
fs 

CO 

O 

< 

CO 



> 

< 



Remark. The cases of existence of the corresponding (tensor) invariants are indicated by gray color in the table. The partial filling corresponds to the uncomplete set of integrals. 
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